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EFFECTS  OF  CRIMPED  FIBER  PATHS  ON  MIXED-MODE  DELAMINATION 
BEHAVIORS  IN  WOVEN  FABRIC  COMPOSITES 


INTRODUCTION 


Woven  fabric-reinforced  polymer  (WFRP)  composites  have  become  an  important 
resource  for  structures  that  require  lightweight  materials  and  are  capable  of  operating  in  and 
surviving  against  severe  dynamic  loading  events  such  as  blast,  ballistic  and  fragment  impacts, 
and  mechanical  shock.  The  applicability  of  WFRP  composites  to  primary  load-carrying 
structures,  however,  is  often  limited  because  of  complexities  involving  localized  mechanics  that 
span  across  multiple  material  length  scales.  These  complexities  influence  the  failure  thresholds 
and  damage  mechanisms  in  WFRP  composites  and  are  difficult  to  simulate  in  numerical  models. 
Specifically,  numerical  simulations  that  employ  homogenization  techniques  at  the  ply  and 
laminate  levels  explicitly  lack  the  critical  details  of  the  fiber  paths  and  the  fiber-matrix  interfaces 
to  address  the  influence  of  woven  fabric  architectures  on  fracture  behaviors  and  interfacial 
delaminations. 

The  research  documented  in  this  report  shows  that  these  behaviors  can  be  investigated 
through  mesoscale  multi-continuum  models,  the  data  from  which  can  contribute  to  the 
identification  of  new  and  advanced  tailoring  methods  to  improve  structural  performance  of 
composites  by  increasing  fracture  resistance,  damage  tolerance,  and  energy  absorption 
capacities. 


PUBLISHED  RESEARCH  RELATED  TO  CHARACTERIZING  FRACTURE 
TOUGHNESS  AND  CRACK  PROPAGATION  PROPERTIES 


This  section  highlights  salient  developments  in  performance  modeling  for  fracture 
toughness  and  crack  propagation  behaviors  of  fiber-reinforced  polymer  (FRP)  composite 
laminates.  A  commonly  used  method  to  perform  numerical  modeling  of  mixed-mode  crack 
growth  in  laminated  composites  is  the  virtual  crack  closure  technique  (VCCT),  originally 
proposed  by  Rybicki  and  Kanninen.1  The  technique  is  based  on  linear  elastic  fracture  mechanics 
(LEFM)  and  the  principle  that  the  energy  released  during  crack  extension  is  the  energy  required 
to  close  the  crack.  Krueger2  provided  in-depth  descriptions  of  VCCT  historical  developments, 
detailed  theoretical  descriptions  for  its  implementation  in  finite-element  analysis  (FEA)  codes, 
and  descriptive  examples  for  multidimensional,  bi-material  interface  applications.  Unlike 
traditional  stress-based  approaches,  which  require  locally  refined  crack-tip  meshes  and 
specialized  quarter-point  crack-tip  elements  to  yield  the  1  /  Vr  singularity  (where  r  is  the  radial 
distance  from  crack  tip),  the  VCCT  provides  an  efficient  alternative  to  numerical  modeling  of 
crack  propagations  by  using  a  kinematic  approach  to  structures  containing  an  initial  crack.  Here 
the  nodal  forces,  opening  and  sliding  displacements,  are  tracked,  and  the  mode-specific  strain 
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energy  release  rates  are  computed  and  compared  to  fracture  toughness  properties  of  the  material 
modeled.  Crack  growth  occurs  when  the  effective  mixed-mode  strain  energy  release  rate 
exceeds  the  material’s  fracture  toughness  in  accordance  with  a  governing  mixed-mode  crack 
growth  criterion.  Additional  considerations  are  described  for  geometric  nonlinearities  and 
mismatched  element  sizes  between  interface  surfaces. 

'J 

Benzeggagh  and  Kenane  conducted  experimental  tests  and  fractography  on  6-mm-thick 
E-glass/epoxy  composite  laminates  subjected  to  mode  I  double-cantilever  beam  (DCB),  mode  II 
end-loaded  split  (ELS),  and  mixed-mode  (I  and  II)  bending  (MMB)  loads  in  displacement 
control.  Benzeggagh  and  Kenane3  presented  a  semi-empirical  criterion  for  describing  the 
mixed-mode  fracture  toughness — known  as  the  BK  mixed-mode  fracture  toughness  criterion. 
Fractography  analysis  using  a  scanning  electron  microscope  identified  the  presence  of  different 
fracture  surfaces  resulting  from  the  mode  I,  mode  II,  and  mixed-mode  loadings.  Mode  I  loading 
produced  cleavage  fractures  in  the  epoxy  and  along  the  epoxy/fiber  interfaces.  Mode  II  loading 
produced  epoxy  fractures  with  shear  hackles  oriented  at  45°  to  the  fiber  paths.  This  orientation 
was  coincident  with  the  plane  of  maximum  principal  stress,  thus  confirming  the  dominance  of  in¬ 
plane  shear  stresses.  The  mixed-mode  loading  produced  fracture  surfaces  with  combined 
cleavage  and  shear  hackle  fractures;  however,  as  the  mode  mixity  increased  with  increasing 
mode  II  contributions,  the  quantity  of  shear  hackles  increased  and  cleavage  fractures  decreased. 

Shindo  et  al.4  investigated  the  mixed-mode  fracture  behaviors  in  woven  glass/epoxy 
laminates  at  several  temperatures:  room  temperature,  liquid  nitrogen  temperature  (77°K),  and 
cryogenic  temperature  using  liquid  helium  (4°K).  Experiments  were  performed  using  the 
ASTM-D-6671-06  MMB  test5;  both  analytical  and  three-dimensional  (3-D)  FEA  approaches 
were  used  to  obtain  strain  energy  release  rates.  The  experimental  results  showed  that 
(1)  unstable  crack  growth  occurred  at  the  77°K  and  4°K  temperatures  and  (2)  the  total  critical 
strain  energy  release  rates  at  room  temperature  and  at  4°K  increased  with  increasing  mixity  from 
mode  II  contributions.  The  FEA  model  treated  the  woven  laminate  as  a  homogeneous, 
orthotropic  material  and  used  the  VCCT  and  BK  mixed-mode  fracture  criterion  to  govern  crack 
propagation.  Homogenization  performed  at  the  ply  level  and  above — a  commonly  used 
approach  for  modeling  WFRP  composites — does  not  allow  the  direct  representation  of  the 
crimped  fiber  paths;  therefore,  the  effects  of  crimp  on  fracture  toughness  cannot  be  simulated. 

Ranatunga6  investigated  delamination  initiation  and  crack  propagation  behaviors  in 
unidirectional  carbon  fiber/epoxy  composite  laminates  using  laminate-level,  homogenized  FEA 
models  that  employed  the  VCCT  and  cohesive  zone  modeling  (CZM)  methods.  The  models 
were  representative  of  the  DCB  test  for  evaluating  mode  I  fracture  toughness.  The  DCB 
specimen  was  similar  to  the  MMB  specimen  and  included  a  film  insert  to  provide  the  initial 
debond  length.  Parameters  influencing  solution  convergence  and  runtimes,  including  mesh 
refinement,  viscous  and  stabilization  energies,  penalty  stiffness  (for  CZM),  and  time-step  size, 
were  investigated.  As  demonstrated,  both  the  VCCT  and  CZM  approach  to  delamination 
initiation  and  propagation  modeling  must  carefully  consider  these  parameters  to  obtain  reliable 
solutions. 

Camanho  and  Davila7  described  the  development  and  use  of  zero-thickness  volumetric 
decohesion  finite  elements  to  simulate  delamination  initiation  and  crack  growth  in  composite 
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laminates.  Their  approach  was  similar  to  that  of  the  CZM.  A  quadratic  softening  law  was 
proposed  for  mixed-mode  loadings.  Both  single-  and  mixed-mode  progressive  delamination 
models  were  developed.  Comparison  of  the  model  results  and  physical  test  data  showed  good 
agreement. 

Motamedi  and  Milani8  discussed  fracture  analysis  methods  for  modeling  delaminations  in 
composite  materials  subjected  to  mixed-mode  loadings.  These  methods  included  the  VCCT, 
CZM,  and  extended  finite-element  methods  (XFEMs)  for  evaluating  nonlinear  fracture  problems. 

Hulton  and  Cavallaro9  investigated  fracture  of  a  thin-walled,  carbon  fiber/epoxy  matrix 
cylinder  subject  to  quasi-static  diametral  compression  using  experiments  and  XFEM  modeling. 
The  cylinder  was  constructed  of  unidirectional  (non-woven)  tapes  oriented  along  the  hoop  and 
longitudinal  axes.  The  XFEM  model  used  ply-level  homogenization.  Crack  initiation  loads  and 
crack  growth  behaviors  were  monitored.  Comparison  of  experimental  and  XFEM  results 
showed  that  the  predicted  primary  fracture  mode  and  crack  initiation  site  matched  those  observed 
from  testing.  Fracture  initiated  via  mode  I  behavior  in  the  outer  longitudinal  ply  propagated 
radially  inward  and  deflected  at  an  inner  90°  ply;  a  mode  I  cohesive  interface  failure  developed 
along  the  hoop  direction.  The  study  demonstrated  the  effectiveness  of  using  XFEM  to  model 
fracture  in  fiber-reinforced  composites  without  requiring  prior  knowledge  of  crack  initiation  and 
propagation  paths. 

Hawkins  and  Haque10  performed  experimental  mode  I  tests  on  10-ply  unidirectional 
carbon/epoxy  laminates  with  and  without  graphene  nanoparticle  reinforcements  (0.1%  by 
weight)  using  the  DCB  specimen  in  displacement-controlled  loading.  The  mode  I  critical  strain 
energy  release  rates  (G/c)  were  measured  and  compared.  It  was  concluded  that  the  three-phase 
carbon/graphene/epoxy  laminate  provided  an  8.3%  increase  in  Gjc  over  the  conventional 
carbon/epoxy  laminates. 

Cavallaro11  experimentally  investigated  the  effects  of  several  fabric  weave  styles  and 
crimp  gradients  on  the  damage  tolerance  and  energy  absorption  capacities  in  20-ply 
Kevlar/epoxy  composite  laminates.  The  experimental  results  demonstrated  that  weave  styles  and 
crimp  gradients  can  be  used  to  alter  the  spatial  and  temporal  behaviors  of  stress  wave 
propagations  and  can  increase  damage  tolerance  and  energy  absorption  capacities  for  static  and 
severe  dynamic  loading  events. 

This  report  documents  the  results  of  research  that  focused  on  the  influence  of  the 
curvilinear  fiber  paths  developed  from  weaving  on  the  fracture  toughness  and  crack  propagation 
behaviors  of  the  fiber/matrix  interfaces  by  performing  two-dimensional  (2-D),  multi-continuum 
FEA  using  VCCT  with  the  BK  mixed-mode  fracture  criterion.  This  research  also  explored  the 
presence  of  any  localized  mixed-fracture  modes  produced  at  these  interfaces  as  a  result  of 
single-mode  global  loadings.  The  presence  of  localized  mixed  modes  can  be  leveraged  as  a 
mechanism  for  increasing  the  fracture  toughness  and  reducing  interfacial  delaminations  in 
WFRP  composites.  Unlike  unidirectional  polymer  composites  in  which  the  fiber  paths  are 
straight,  WFRP  composites  have  increased  damage  tolerance,  possibly  resulting  from  localized 
mixed-mode  fracture  behaviors.  The  findings  of  this  research  will  provide  new  and  advanced 
methods  for  improving  the  structural  performance  of  WFRP  composites. 
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MODEL  DEVELOPMENT 


Multi-continuum  mesoscale  models  were  developed  using  ABAQUS12  FEA  software  to 
quantify  the  primary  and  secondary  strain  energy  release  rates  that  are  present  during  pure 
mode  I,  mode  II,  and  mixed-mode  global  loadings  on  delaminations  of  the  fiber/matrix 
interfaces.  Figure  1  illustrates  the  conventional  fracture  modes  considered  for  this  modeling 
effort.  The  component  strain  energy  release  rates  were  obtained  using  VCCT;  the  mixed-mode 
responses  were  resolved  using  the  BK  fracture  criteria.  The  variations  of  localized  mode  I  and 
mode  II  strain  energy  release  rates  were  captured  over  a  one-half  crimp  wavelength  and 
compared.  The  numerical  results  provide  confirmation  that  fracture  behaviors  are  dependent  on 
the  crimped  fiber  paths  and  that  mixed-mode  fracture  behaviors  exist  at  the  mesoscale  for  single- 
and  mixed-mode  global  loadings. 


Mode-1  Mode-ll  Mode-Ill 

(Opening)  (In-Plane  Shear)  (Out-of-Plane  Shear) 


Figure  1.  Fracture  Modes 


Three  laminate  models  were  developed  at  the  mesoscale  with  each  laminate  constructed  of 
three  plain-woven  plies.  The  plies  for  each  individual  laminate  were  consistent;  however,  the 
three  models  contained  different  crimped  fiber  paths  and  were  identified  in  terms  of  their  warp 
and  weft  crimp  distributions  as  fully  balanced  (FB),  semi-balanced  (SB)  and  unbalanced  (UB)  as 
shown  in  figure  2.  The  laminate  models  were  loaded  in  a  manner  representative  of  the  MMB  test 
fixture. 


3-Ply  Plain-Woven  Laminate  FEA  Models 

Fully-Balanced  (FB)  Crimp 


Semi-Balanced  (SB)  Crimp 


Weft  Fibers 


Warp  Fibers 


Unbalanced  (UB)  Crimp 


Figure  2.  Three  Laminate  Models  Constructed  of  Three-Ply  Plain-Woven  Fabrics 

with  Different  Crimped  Fiber  Paths 
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The  fiber  and  matrix  materials  were  idealized  as  linearly  elastic  and  isotropic  using  the 
properties  listed  in  table  1 .  A  more  definitive  representation  of  actual  fibers  modeled  as 
homogenized  yam  forms  would  suggest  using  a  transversely  isotropic  behavior;  however, 
sufficient  material  property  data  were  not  available  to  describe  the  fiber  yam  forms.  Cross 
sections  of  the  weft  fiber  yams  were  modeled  as  ellipsoidal.  The  separation  distances  between 
the  warp  and  weft  fibers  were  consistent  across  all  three  laminate  models  (see  table  2  for 
additional  model  geometry  details). 


Table  1.  Material  Properties 


Material 

E  (psi) 

V 

G  (psi) 

Fiber 

10.0  x  106 

0.33 

3.74  x  10b 

Matrix 

1.0  x  106 

0.27 

3.94  x  105 

Table  2.  Model  Geometry  Details 


Parameter 

Value 

Model  Length 

1.2512  in. 

FB  Laminate  Thickness 

0.1080  in. 

SB  Laminate  Thickness 

0.1347  in. 

UB  Laminate  Thickness 

0.1620  in. 

Warp  Thickness 

0.0126  in. 

Maximum  Weft  Thickness 

0.0126  in. 

Maximum  Weft  Width 

0.0548  in. 

Weft  Width/Thickness  Ratio 

4.3492 

Weft  Density 

6.394  yams/in. 

One-Half  Crimp  Wavelength  (A/2) 

0.0780  in. 

Initial  Debond  Length  ( ao ) 

0.1560  in. 

The  MMB  test  fixture,  initially  designed  by  Reeder,13  was  represented  in  the  model 
shown  in  figure  2  to  apply  the  global  mode  I,  global  mode  II,  and  global  mixed-mode  loadings 
(modes  I  and  II).  Note  that  this  fixture  and  the  2-D  models  do  not  support  mode  III  loadings. 

The  degree  of  mode  mixity  was  controlled  by  varying  the  location  of  the  applied  force  as 
shown  in  figure  3.  The  force  was  applied  to  the  loading  bar  at  a  constant  rate  with  respect  to 
time.  This  type  of  loading  is  referred  to  as  “force-control  loading”  and  is  an  alternative  to 
displacement-controlled  loading,  in  which  the  displacement  is  applied  at  a  prescribed  rate  with 
respect  to  time.  Force-control  loading  was  the  preferred  type  for  this  research  because  it  is  most 
representative  of  loading  events  on  realistic  structures. 
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•  Force  (-G- )  is  applied  individually  at  constant  rate  over  time 

•  When  force  is  applied  at  X  =L,  global  mode-1  fracture  results. 

•  When  force  is  applied  at  X  =3L/4,  global  mixed  mode  (1,11)  fracture  results. 

•  When  force  is  applied  at  X  =L/2,  global  mode-11  fracture  results. 


Figure  3.  Method  of  Loading  and  Boundary  Conditions  for  Mixed-Mode  Bending  Test 
Simulations  of  Three-Ply  Woven-Fabric  Composite  Laminate 

Plane  stress  (2-D)  continuum  elements  of  type  CPS4  (four-node  plane  stress  element) 
were  used  to  represent  the  warp  and  weft  fibers  and  matrix  components.  A  plane  stress 
formulation,  in  contrast  to  plane  strain,  was  chosen  because  the  model  considers  the  shallow 
depth  of  a  single  warp  yam  with  respect  to  the  laminate  depth  (refer  to  global  Z-axis  in  figure  3). 
Contact  surfaces  were  defined  along  the  fiber/matrix  interfaces  for  each  ply  and  were  assigned 
tied  kinematic  restraints.  Note  that  the  initial  debond  surfaces  for  the  woven  composites  were 
not  straight  lines  (or  flat  planes  in  3-D)  as  those  used  for  unidirectional  composites;  rather,  they 
were  curved  surfaces  along  the  central  ply’s  bottom  warp  (longitudinal)  fiber/matrix  interface. 

No  tied  contact  restraints  were  assigned  to  nodes  along  the  initial  debond  length,  as 
shown  in  figure  4.  The  absence  of  tied  contact  restraints  at  nodes  above  and  below  the  initial 
debond  regions  enabled  the  unbonded  surface  region  nodes  to  move  independently.  The  tied 
contact  restraints  ahead  of  the  crack  tip  were  eliminated  when  the  total  strain  energy  release  rate 
Gtotai  reached  the  effective  critical  strain  energy  release  rate  Gc_eff.  Once  this  occurred,  strain 
energy  was  released  and  crack  growth  developed.  Failure  was  based  on  LEFM  through  the  use 
of  VCCT.  The  VCCT  supports  crack  propagation,  but  it  does  not  predict  crack  initiation.  An 
initial  crack  must  be  present.  The  mode  I  and  mode  II  strain  energy  release  rates  that  determine 
crack  propagation  are  reported  when  using  VCCT.  Another  commonly  used  modeling  approach 
to  fracture  analysis  of  composite  laminates  is  CZM  using  a  traction-separation  law  (TSL).  The 
CZM  incorporates  cohesive  surfaces  with  a  TSL  that  defines  a  combined  strength/damage 
initiation,  damage  evolution  criteria,  and  penalty  stiffness  to  model  both  crack  initiation  and 
crack  propagation;  however,  it  does  not  report  the  individual  strain  energy  release  rates.  This 
research  therefore  used  the  VCCT  with  a  representative  MMB  test  specimen  and  loading. 


Desired  minimum  crack 
growth  length  =  M2 


Figure  4.  Initial  Debond  Length  and  Desired  Minimum  Crack  Growth  Size 

(Red  path  represents  lower  fiber/matrix  interface  of  the  central  warp  fiber 
where  crack  growth  is  permitted.) 
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Displacement  boundary  conditions  included  a  pinned  support  at  the  laminate’s  lower  left 
end  and  roller  support  at  the  laminate’s  lower  right  end.  A  rigid  loading  bar  contacted  the  upper 
surface  of  the  laminate  at  the  midspan  (surface-to-surface  contact)  and  at  the  right  end  (tied 
contact).  The  effects  of  geometric  nonlinearities  from  large  displacements  and  contact  between 
the  loading  bar  and  the  upper  surfaces  of  the  laminates  were  included. 

The  BK  mixed-mode  fracture  toughness  criterion  was  developed  by  Benzeggagh  and 
Kenane3  and  is  described  by  equation  (1)  for  2-D  cases.  For  3-D  cases,  the  BK  criterion  includes 
a  mode  III  strain  energy  release  rate  term  as  shown  in  equation  (2). 


G C_eff  —  Gic  +  (GIIC  GIC) 


total 


,  and 


(1) 


Gc_eff  —  Gjc  +  (GIIC  GIC ) 


G„  +  Gm 


J  total 


(2) 


where 

Gi ,  Gu ,  Gm  are  the  modes  I,  II,  and  III  strain  energy  release  rates,  respectively, 

Gic ,  Guc ,  Guic  are  the  modes  I,  II,  III  fracture  toughness,  respectively, 

G total  =Gi+  Gu  +  Gm, 

77  is  a  curve  fit  parameter  that  was  shown  to  vary  with  the  matrix  ductility, 

Gc_eff  is  the  mixed-mode  fracture  toughness,  and 

strain  energy  release  rates  and  fracture  toughness  terms  are  expressed  in  (in.-lb/in.  )  units. 

The  open  literature  contains  minimal  data  regarding  the  mode  I  and  mode  II  critical  strain 
energy  release  rates  for  laminated  woven  composites.  The  mode  I  and  mode  II  critical  strain 
energy  release  rates  and  material  parameter  77  reported  by  Ranatunga6  for  unidirectional  AS4 
carbon  fiber/350 1-6  epoxy  laminates  were  used  in  this  research  after  the  units  of  measure  were 
converted  to  the  English  system  of  units. 


COMPARISON  OF  INTERFACIAL  FRACTURE  BEHAVIORS 


Solutions  were  obtained  for  each  laminate  type  at  each  of  the  three  applied  force 

f-  f" 

locations  {X  =  L,X  =  ILIA,  and  X  =  LI 2)  (see  figure  3).  Figure  5  shows  the  effects  of  loading 
position  on  the  global  Sxx  stress  component  for  the  SB  laminate.  The  applied  loads  were 
increased  until  interfacial  crack  growth  occurred  for  a  minimum  of  distance  of  XU.  Note  that  the 
cleavage  deformations  shown  in  the  top  image  of  figure  5  correspond  to  the  mode  I  fracture.  The 
shear-dominated  deformations  are  shown  in  the  inset  view  for  the  mode  II  fracture  case. 
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Figure  5.  Contour  Plots  of  Sxxfor  the  SB  Laminate  for  Applied  Load  Locations 

atX=L,X  =  3L/4,  X  =  L/2 

(Note  the  extensive  crack  growth  for  load  applied  at  X  =  L  (mode  I  global  loading)). 


The  contour  plots  in  figures  6-8  illustrate  the  values  of  Gtotai/Gc_eff,  Gi,  and  Gn  for  the  FB 
laminate  with  loading  applied  at  the  X  =L,  X=3L/4,  and  X  =L/2  locations.  The  strain  energy 
release  rates  Gy ,  Gn ,  Gtoto/,  the  Glolai/Gc_eff  ratios,  and  fracture  toughnesses  Gc_eff  for  the  FB,  SB, 
and  UB  laminates  were  plotted  as  functions  of  position  along  the  fiber  paths,  as  shown  in 
figures  9-11.  The  FB  laminate  exhibited  the  greatest  variations  in  strain  energy  release  rates  and 
fracture  toughness  values  along  the  fiber  paths  followed  by  the  SB  and  UB  laminates, 
respectively.  Most  notable  were  the  regions  of  zero  G/  contributions  to  the  mixed-mode  fracture 
toughness  Gc  eff  in  the  FB  and  SB  laminates  occurring  for  loads  applied  at  the  X  =  LI 2  (global 
mode  II)  and  X  =  3T/4  (global  mixed  modes  I  and  II)  locations.  These  regions  were  coincident 
with  large  fiber-path  angles  and  experienced  localized  conversions  of  G/  to  Gn  strain  energy 
release  rates  and  both  Gtotai  and  Gc_eff  increased  by  nearly  3:1  with  respect  to  their  values  at  the 
X=0  fiber-path  positions. 

Reviewing  the  Gc_eff  curves,  it  can  be  seen  that  delaminations  will  develop  earlier  in  the 
FB  and  SB  laminates  than  in  the  UB  laminate;  however,  the  FTCM  arrested  further  crack  growth 
until  Gtotai  exceeded  the  higher  Gc_eff  values  resident  at  the  X  =  0.025  through  X  =  0.045  fiber- 
path  positons. 
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Figure  6.  Contour  Plots  of  Gtota/Gc_eff,  Git  and  Gnfor  the  FB  Laminate  at 
Warp  Fiber/Matrix  Interface,  Load  Applied  at  X  =  L 
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Figure  7.  Contour  Plots  of  Gtota/Gc  eff,  GIy  and  Gnfor  the  FB  Laminate  at 
Warp  Fiber/Matrix  Interface ,  Load  Applied  at  X  =  3L/4 
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Figure  8.  Contour  Plots  of  Gtota/Gceff>  6/,  and  Gnfor  the  FB  Laminate  at 
Warp  Fiber/Matrix  Interface ,  Load  Applied  at  X  =  L/2 


X-Distance  From  Initial  Crack  Tip  (in) 

(a)  Ratios  of  GtotJGc_eff  Versus  X-Distance 


FULLY  BALANCED 


X-Distance  From  Initial  Crack  Tip  (in) 

(b)  Strain  Energy  Release  Rates  Versus  X-Distance 


FULLY  BALANCED 


X-Distance  From  Initial  Crack  Tip  (in) 

(c)  Gc _eff  Versus  X-Distance 


Figure  9.  Plots  of  Gtota/Gc_eff  Ratio,  Strain  Energy  Release  Rates  Gi,  Gn,  Gtotai,  and  Fracture 
Toughness  Gc  eff  for  the  FB  Laminate  Versus  X-Distance  from  Initial  Crack  Tip 
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SEMI  BALANCED 

1.4  r  (G/c=0.457,  G//c=3.141,  G„,c=0.0,  17=1.80) 


(Shown  at  Final  Time  Step) 


0.00  0.01  0.02  0.03  0.04  0.05 


0.07  0.08  0.09 


X-Distance  From  Initial  Crack  Tip  (in) 

(a)  Ratios  of  GtotailGc_eff  Versus  X-Distance 


4.0 


SEMI  BALANCED 

(G/C=0.4S7,  G(/c=3.141,  Gmc=0.0,  7=1.80) 
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(b)  Strain  Energy  Release  Rates  Versus  X-Distance 


SEMI  BALANCED 


(G/c=0.457,  G;;c=3.141, 


G///c=0.0,  7=1.80) 
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Mode-I  Global  Loading 
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(c)  Gc_  eff  Versus  X-Distance 


Figure  10.  Plots  of  Gtota/Gc_eff  Ratio,  Strain  Energy  Release  Rates  Gi,  Gn,  Gtotai,  and  Fracture 
Toughness  Gceff  for  the  SB  Laminate  Versus  X-Distance  from  Initial  Crack  Tip 


13 


(a)  Ratios  of  GtotJ Gc_eff  Versus  X-Distance 


UNBALANCED 

(G^O.457,  G,k=  3.141,  G,Wt= 0.0,  <7=1.80) 


— G total  AT  L  -»-GI  AT  L  — Gil  AT  L 

— Gtotal  AT  3L/4  -»-GI  AT  3L/4  — Gil  AT  3L/4 

- Gtotal  AT  1/2  -*-GI  AT  L/2  -»-GII  AT  L/2 


Xt  tTi 


X-Distance  From  Initial  Crack  Tip  (in) 


(b)  Strain  Energy  Release  Rates  Versus  X-Distance 
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(c)  Gc_e#  Versus  X-Distance 


Figure  11.  Plots  of  Strain  Energy  Release  Rates  Gi,  Gn,  Gtotai,  and  Fracture  Toughness 
Gc_eff  for  the  UB  Laminate  Versus  X-Distance  from  Initial  Crack  Tip 
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Table  3  lists  the  maximum  predicted  values  of  Gc_eff  along  the  fiber  paths  and  the  values 
corresponding  to  the  X=  0  fiber  path  position  for  each  laminate  and  global  loading  mode.  The 
Gc_eff  values  were  also  plotted  in  figure  12. 


Table  3.  Effective  Critical  Strain  Energy  Release  Rates,  Gc _efffor  0  <x  < A/ 2 


Global  Loading  Mode 

Laminate  Type 

Gc  eff  at  X  —  0 
(in.-lb)/in.2 

Maximum  Gc  eff 
(in.-lb)/in.2 

I 

FB 

0.859 

1.770 

SB 

0.501 

0.555 

UB 

0.505 

0.870 

II 

FB 

0.881 

3.141 

SB 

1.505 

2.907 

UB 

2.190 

2.640 

Mixed  Modes(I,  II ) 

FB 

1.020 

3.141 

SB 

1.680 

2.990 

UB 

2.160 

2.760 

Global  Mode-ll  Loading  (x=  L/2) 


3.5 


At  X=0  ■  Max 


Figure  12.  Plots  of  the  Effective  Critical  Strain  Energy  Release  Rates,  Gceff 
Versus  Laminate  and  Global  Loading  Type 


15 


The  observed  strain  energy  release  rate  conversion  mechanism  was  further  investigated 
by  examination  of  the  stress  components  along  the  fiber  paths.  A  stress  transformation  was 
performed  using  equation  (3)  to  resolve  the  global  coordinate  system  to  a  local  coordinate 
system.  The  local  system  provided  the  interlaminar  tangential,  normal,  and  shear  stress 
components  Sue,  Sly,  and  Sucy,  respectively,  along  the  fiber/matrix  interfaces.  Figures  13-15 
show  the  global  and  interlaminar  stress  components  at  the  immediate  time  step  preceding  crack 
growth.  Note  that  the  global  and  local  coordinate  systems  for  the  UB  laminate  model  were 
identical;  therefore,  the  interlaminar  stresses  were  equal  to  the  global  stress  components. 
Interfacial  delaminations  will  initiate  when  the  interlaminar  normal  stress  Sly  is  tensile  in  sense 
and/or  the  interlaminar  shear  stress  Slxy  exceed  their  respective  interface  strength  levels. 


Slx  1 

cos2  6 

sin2Q 

—2cos6sin0 ' 

(  Sx 

II 

to1 

sin2d 

cos26 

IcosOsinO 

WIT 

Slxy) 

-cosdsinO 

—cosdsinO 

cos26—sin26- 

where  0  is  measured  from  the  tangent  vector  along  the  fiber  path  to  the  global  x-axis  and  is 
expressed  in  radians,  and  Sx,  Sy,  Sxy,  Slx,  Sly, and  Slxy  are  expressed  in  units  of  psi. 
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FULLY  BALANCED,  LOAD  AT  L 
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X-Distance  From  Initial  Crack  Tip  (in) 

<— 

(a)  Load  Applied  at  X  =  L 


(b)  Load  Applied  at  X  =  3Z/4 


G,c=0.457,  G//c=3.141,  G„,c=0.0,  77  =1.80 

(STRESSES  SHOWN  PRIOR  TO  CRACK  GROWTH) 


J_ I_ I_ I_ I_ I_ L 


X-Distance  From  Initial  Crack  Tip  (in) 

(c)  Load  Applied  at  X  =  L/2 

Figure  13.  Plots  of  Global  and  Local  Stress  Components  Versus  X-Distance  from  Initial 
Crack  Tip  Shown  Immediately  Preceding  Crack  Growth  for  the  FB  Laminate 


17 


16,000 
14,000 
12,000 
10,000 
^  8,000 
£  6,000 
^  4,000 
2,000 
0 

-2,000 


SEMI  BALANCED,  LOAD  AT  L 

G/c=0.457,  G„c=3.141,  G,„c=0.0,  77=1.80 
(STRESSES  SHOWN  PRIOR  TO  CRACK  GROWTH) 
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(a)  Load  Applied  at  X  =  L 


0.09 


SEMI  BALANCED,  LOAD  AT  3L/4 


(b)  Load  Applied  at  X  =  3L/4 


SEMI  BALANCED,  LOAD  AT  L/2 


(c)  Load  Applied  at  X  =  LI 2 


Figure  14.  Plots  of  Global  and  Local  Stress  Components  Versus  X-Distance  from  Initial 
Crack  Tip  Shown  Immediately  Preceding  Crack  Growth  for  the  SB  Laminate 
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UNBALANCED,  LOAD  AT  L 
G/c=0.457,  G„c=3.141,  G/;/c=0.0,  17=1.80 
(STRESSES  SHOWN  PRIOR  TO  CRACK  GROWTH) 
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<— 

(a)  Load  Applied  at  X  =  L 


UNBALANCED,  LOAD  AT  3L/4 
G,c=0.457,  G</c=3.141,  G,„c=0.0,  #7=1.80 
(STRESSES  SHOWN  PRIOR  TO  CRACK  GROWTH) 
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(b)  Load  Applied  at  X  =  3Z/4 


UNBALANCED,  LOAD  AT  L/2 
G,c=0.457,  G„c=3. 141,  G„/c= 0.0,  7=1.80 

(STRESSES  SHOWN  PRIOR  TO  CRACK  GROWTH) 
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(c)  Load  Applied  at  X  =  L/2 


Figure  15.  Plots  of  Global  Stress  Components  Versus  X-Distance  from  Initial  Crack  Tip 
Shown  Immediately  Preceding  Crack  Growth  for  the  UB  laminate 
(Note  that  the  global  and  local  stress  components  were  the  same  for  the  UB  laminate.) 
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The  applied  loads  immediately  prior  to  crack  growth  and  the  peak  local  stresses  in  the 
vicinity  of  the  initial  crack  tip  (atX =  0.004  inch)  immediately  prior  to  crack  growth  are  listed  in 
table  4.  The  loads  were  observed  to  increase  from  the  mode  I  to  mode  II  to  the  mixed-mode  case 
as  expected  because  of  the  differences  in  the  GIc  and  Gnc  values  used.  In  a  similar  fashion,  the 
loads  were  also  observed  to  increase  with  laminate  type  in  increasing  order  from  the  FB  laminate 
to  the  SB  laminate  to  the  UB  laminate.  Recall,  however,  that  the  laminate  thickness  for  each  type 
increases  accordingly  as  shown  in  table  2.  The  normal  and  shear  stress  components  generally 
dominated  the  three  types  of  global  loadings.  (Note  that  the  local  and  global  coordinate  systems 
were  identical  for  the  UB  laminate  and,  as  such,  the  corresponding  stress  values  in  table  4  were 
simply  designated  using  the  local  stress  nomenclature.)  The  shear  strain  energies  served  to 
increase  Gc _ejf  for  the  FB  and  SB  laminates  when  they  were  subjected  to  mixed-mode  and  mode  II 
global  loadings.  This  phenomenon  was  also  observed  by  Shindo.4  The  local  tangential  stress 
components  Slx  were  the  smallest  values  for  the  mode  II  and  mixed-mode  loadings  across  each 
laminate  style.  The  maximum  interface  stresses  achieved  just  prior  to  crack  growth  were 
generally  observed  for  the  UB  laminate  for  each  of  the  three  loading  modes.  For  this  laminate 
and  the  mode  II  and  mixed-mode  loadings,  the  shear  stresses  were  the  maximum  components 
developed,  followed  by  the  normal  and  tangential  components,  respectively.  The  maximum  shear 

f-  f- 

stresses  achieved  for  the  mixed  mode  ( X  =  321/4 )  and  mode  II  (X  =  LI 2)  loadings  of  the  UB 
laminate  were  nearly  identical  as  the  mode  I  strain  energies  developed  in  the  mixed  mode  case 
had  minimal  influence  on  the  stress  distributions  near  the  crack  front. 

The  stress  components  in  table  4  are  especially  useful  for  defining  TSLs  for  modeling 
woven  composite  laminates  with  CZMs;  however,  it  is  important  that  the  TSLs  are  defined  using 
the  local  stress  components  and  strain  energy  release  rates  to  ensure  proper  calculations  of 
interfacial  delaminations. 
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Table  4.  Applied  Loads  and  Peak  Local  Stress  Components 
(Near  Crack  Tip  at X=  0.004  inch)  Prior  to  Crack  Growth 


FB  Laminate 


Mode  I 

Mode  II 

Mixed  Modes 

Load  (lb) 

53.46 

129.00 

204.60 

Slx 

(psi) 

3223.92 

2717.98 

2916.93 

Sly 

(psi) 

6484.23 

7618.80 

7600.95 

Slxy 

(psi) 

4124.33 

6425.04 

7478.75 

SB  Laminate 


Mode  I 

Mode  II 

Mixed  Modes 

Load  (lb) 

62.75 

233.6 

381.5 

Slx 

(psi) 

3438.53 

2,359.75 

2124.83 

Sly 

(psi) 

6003.55 

6789.23 

6716.27 

Slxy 

(psi) 

3718.03 

10898.10 

11384.30 

UB  Laminate 


Mode  I 

Mode  II 

Mixed  Modes 

Load  (lb) 

78.52 

325.10 

467.00 

Slx 

(psi) 

4770.03 

4775.18 

4842.57 

Sly 

(psi) 

7586.65 

7915.91 

8103.26 

Slxy 

(psi) 

4486.59 

14890.40 

14917.10 

The  VCCT  also  reports  the  nodal  debond  times  (DBTs)  corresponding  to  the  nodes  along 
the  slave  cohesive  surfaces.  Crack  growth  stability  was  examined  by  plotting  DBT  for  each 
laminate  and  loading  type  versus  position  along  the  fiber  paths  as  shown  in  figure  16.  Unstable 
crack  growth  occurred  when  the  slopes  of  the  DBT  versus  position  approach  zero.  The  FB 
laminate  demonstrated  the  highest  levels  of  crack  growth  stability  followed  by  the  SB  and  UB 
laminates,  respectively;  however,  the  FB  and  SB  laminates  showed  varying  levels  of  crack 
growth  stability  that  were  dependent  on  position  along  the  crimped  fiber  paths.  The  highest 
crack  growth  stability  levels  for  these  two  laminates  were  observed  near  regions  of  maximum 
change  in  slope  along  the  fiber  paths  (approximately  0.025  inch);  however,  crack  growth  became 
completely  unstable  once  crack  lengths  exceeded  0.040  inch.  The  UB  laminate,  which  initially 
exhibited  stable  crack  growth  up  to  a  crack  length  equal  to  1.2e-03  inches  became  completely 
unstable  at  larger  crack  lengths. 
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FULLY  BALANCED 


SEMI  BALANCED 


UNBALANCED 


Figure  16.  Plots  of  Nodal  Debond  Times  for  the  FB  (a),  SB  (b),  and  UB  (c)  Laminates 

Versus  X-Distance  from  Initial  Crack  Tip 


SUMMARY  AND  CONCLUSIONS 


Numerical  models  were  conducted  at  the  mesoscale  for  three  variations  of  a  plain-woven 
fabric  composite  subjected  to  mode  I,  mode  II,  and  mixed-mode  global  loadings.  The  effects  of 
crimped  fiber  paths  on  the  effective  critical  strain  energy  release  rates  were  investigated  using  the 
virtual  crack  closure  technique  with  the  Benzeggagh-Kenane  mixed-mode  fracture  criterion. 
Results  demonstrated  there  were  considerable  effects  for  the  FB  and  SB  laminates  when  they 
were  subjected  to  mixed-mode  and  mode  II  global  loadings.  Increasing  contributions  from  the 
localized  mode  II  strain  energy  release  rates  along  the  crimped  fiber  paths  served  to  increase  the 
mixed-mode  fracture  toughness  Gc_e/f,  as  also  reported  by  Shindo.4 

In  addition  to  the  fiber  bridging,  the  observed  fracture  toughness  conversion  mechanism 
resulting  from  the  curvilinear  fiber  paths  with  GUc  >  GIc  provides  further  insight  into  how  woven 
composites  exhibit  greater  damage  tolerance  thresholds  than  do  unidirectional  composites.  Under 
force-control  loading,  the  FB  laminate  exhibited  the  highest  crack  growth  stability,  followed  by 
the  SB  and  UB  laminates;  however,  the  UB  laminate  in  general  supported  the  highest  stresses 
near  the  vicinity  of  the  initial  debond  region  prior  to  crack  propagation,  provided  virtually  no 
crack  growth  stability. 

Results  of  this  investigation  demonstrate  that  increased  delamination  resistance  and  crack 
growth  stability  can  be  achieved  through  careful  selection  of  the  woven  fabric  architectures  by 
further  consideration  of  the  mixed-mode  fracture  behaviors  at  the  fiber/matrix  interfaces. 

Future  studies  will  conduct  experiments  to  provide  measured  properties  for  candidate  fiber 
materials  as  their  constitutive  behaviors  are  expected  to  be  transversely  isotropic. 
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